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■^ . Abstract: Possibly the only unambiguous verification that gamma-ray bursts (GRBs) 

p ! are at cosmological distances would be the observation of multiple images of a gravita- 

O 



tionally lensed burst. Each images would arrive at a different time, but exhibit identical 
H I light curves. We improve upon previous calculations of GRB lensing by using revised cos- 

^ ■ mological burst models based upon better knowledge of burst spectra, and we consider 

several sets of cosmological parameters. Only lensing by the known population of galaxies 
is considered, and we compute the lensing rate and amplification bias using three differ- 
r> I ent lens models. Assuming that 800 bursts per year are above the detection threshold of 

c^ ' the BATSE experiment on the Compton Gamma-Ray Observatory, we predict one lens- 

ing event every 1.5-25 years, with a median time delay between images of about 7 days. 
However, BATSE has an instrumental efficiency of about 12% for detecting double bursts, 
decreasing the detection rate to fewer than one per century for some cosmological parame- 
ters. This inefficiency can be overcome partially if cases of quadruple imaging by elliptical 
potentials are common, but simulations show that lensing by elliptical potentials enhances 
the detection rate by only about 20%. We conclude that it is unlikely that the BATSE 
experiment will detect a lensed burst, but that a future experiment, which is either much 
more sensitive or has a much more efficient duty cycle, could. 
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1. INTRODUCTION 

Data collected by the Burst and Transient Source Experiment (BATSE) on the Compton 
Gamma-Ray Observatory (CGRO) have generated interest in the cosmological origin sce- 
nario for gamma-ray bursts (GRBs) (see Blaes 1994 for a review of various scenarios). In 
particular, the isotropic distribution of bursts is not compatible with any known galactic 
distribution (disk or halo) of objects (Mao & Paczyhski 1992b). The GRBs must originate 
beyond 18kpc from the galactic center or be part of the solar system {e.g., in the Oort 
cloud) to have the required isotropy. Unfortunately, the flux distribution of bursts alone 
cannot distinguish among the local, extended halo, and cosmological possibilities (Lubin 
& Wijers 1993). 

Since BATSE's spatial resolution (~ 4°) is insufficient to unambiguously associate 
individual bursts with individual high redshift galaxies, perhaps the only definitive proof of 
the cosmological origin of GRBs that BATSE could provide is a detection of a gravitational 
lensing event. A lensed GRB would show multiple images, but unlike quasars, they would 
not be identified as a spatially resolved pair (or quadruplet) . Rather, they would appear as 
two bursts that arrive at different times and are located within the same positional error 
box on the sky. They will exhibit the same light curve, but probably will have different 
amplitudes. 

Since the first suggestion that GRBs might be lensed (Paczyhski 1986), several authors 
have presented detailed calculations of GRB lensing. Mao (1992b) considered lensing by 
the known population of galaxies. He found that GRBs should be lensed at a rate of 
about one per thousand, which is comparable to the rate of quasar lensing (Ostriker & 
Vietri 1986). Other authors have considered lensing by a significant cosmological density of 
point masses (Mao 1992b; Blaes and Webster 1992). Computations of GRB microlensing 
(Paczyhski 1987; Mao 1993) and even "femtolensing" (Gould 1992) have been considered. 
Narayan and Wallington (1992a) focused on what can be learned about the lens itself if a 
lensed burst is observed. Nemiroff et al. (1993) have examined the first 44 GRBs detected 
by BATSE for evidence of lensing. The absence of any lensed bursts in this sample is used 
to set rather loose constraints upon the density of lO^-lO^M© objects in the universe. 

Considering the predicted rate of GRB lensing, it is not surprising that no lensing 
events were found in the small sample examined by Nemiroff et al. (1993). Indeed, if 
only galaxies act as lenses, none should be present in the sample of more than 300 bursts 
catalogued so far. Furthermore, even if there were a lensing event in the catalogued data, 
Wambsganss (1993) has claimed that the amount of noise present in the BATSE light 
curves may overwhelm any hope of identifying it. Nowak and Grossman (1994) find that 
this is excessively pessimistic. They conclude that a large fraction of bursts that are bright 
enough to trigger the BATSE detector (5.5cr above the background noise) are also bright 
enough to have light curves that can be distinguished from other distinct bursts. Since 
previous calculations suggest that in a few more years of data there should be one lensing 
event if GRBs are cosmological, searching for lensed bursts may be an exercise with a 
potentially large return. 

The purpose of this paper is to revise and expand upon the computations of lensing 
by the known population of galaxies. In particular, we find a significantly lower lensing 
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rate than did Mao (1992b), for two important reasons. First, he used a cosmological 
burst model where bursts can be seen to redshifts farther than are hkely. According to 
recent models using improved burst spectra (Fenimore et al. 1993; Tamblyn & Melia 1993, 
hereafter TM93), bursts are not likely to be seen beyond ^ ~ 1, whereas Mao (1992b) used 
a model where they can be seen to 2 ~ 1.5 (Mao and Paczyhski 1992a). Second, Mao did 
not use the appropriate selection criterion for the occurrence of a lensing event, namely 
whether or not the second brightest burst is observable. In §2 we construct cosmological 
models of GRBs using the improved spectrum of TM93 for cosmologies (O, A) — (1.0,0), 
(0.1,0), (0.1,0.9). In §3 we use these three cosmological models to compute the rate of 
observable double bursts, using the population of known galaxies, modeled as singular 
isothermal spheres, as lenses. In addition we consider the lensing rate if the lenses are 
non-singular isothermal spheres or have finite mass and no dark, isothermal halos. The 
lensing rate is high enough that a lensing event has a significant chance of occurring during 
BATSE's lifetime. We discuss how BATSE's instrumental inefficiency for detecting lensing 
events defeats this moderately optimistic rate. In §4 we consider lensing by an elliptical 
potential, and derive from simulations the probability for multiple imaging with two, three, 
or four observable components. If four component bursts are sufficiently common, then 
the instrumental inefficiency can be partly overcome. We demonstrate that lens ellipticity 
enhances the detection rate of lensing events by a small amount. We summarize and discuss 
the implications of our results for the BATSE experiment and potential future experiments 
in 55. 



2. THE GRB POPULATION 

Analysis of the flux distribution of GRBs (Meegan et al. 1992) shows that the observed 
(y/Fmax) statistic (Schmidt 1968) is not compatible with a homogeneous distribution. 
(This statement presupposes that the flux decreases as the inverse square of the distance.) 
Although the (V^/Knax) statistic can be reconciled with a local population that is absent at 
great distances, the observed GRB isotropy is inconsistent with a local galactic halo pop- 
ulation (unless the core radius of the distribution is anomalously large, Mao & Paczyhski 
1992b) and with a local population associated with the Oort cloud, Maoz 1993, Clarke 
et al. 1994). Another possibility is that the inverse square law does not hold, as would be 
the case if bursts are at cosmological distances {z > 0.3) such that the curvature of the 
universe is important. Several authors (Dermer 1992; Mao & Paczyhski 1992a; Piran 1992; 
Wickramasinghe et al. 1993; Fenimore et al. 1993; TM93) have shown that such departures 
from Euclidean geometry can account for the GRB flux distribution. The observation that 
fainter bursts have longer duration (Norris et al. 1994) is consistent with cosmological time 
dilation of distant bursts, and supports the cosmological scenario. 

We assume GRBs originate at cosmological distances, and that the population can 
be described by standard candle luminosity Lq and a constant burst rate per comoving 
volume uq. The maximum distance to which a burst can be seen measures the departure 
from Euclidean geometry, and thus, Lq is chosen to reproduce the observed (F/Fmax)- 
Based on 207 bursts in the BATSE burst catalog (Fishman et al. 1992) up to March 5, 
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1992 (overwrites are excluded), we compute (F/Fmax) = 0.32, in agreement with that from 
the shghtly smaUer sample in Meegan et al. (1992). The burst rate per unit volume no is 
chosen to reproduce the observed burst rate, A^ ~ 800 yr~^. 

The details of this calculation depend on the GRB spectrum. We adopt the double 
power-law model of TM93, which is based on the analysis of burst spectra by Band et al. 
(1993). The photon spectrum [i.e., photons per unit time per unit energy) is described by 

Se{E) = I ^^ai-a'a^aa-l^ ^ > ^^^ (1) 

where Ei, = 300 keV is the break in the power-law slope. The power law exponents are 
cti = 1, a2 = —1. The coefficient A is calibrated by equating the energy emitted in the 
bandpass 50-300 keV (the bandpass in which BATSE bursts are observed) to the restframe 
luminosity, Lq. 

The photon flux observed in the 50-300 keV bandpass depends on the distance to the 
burst, measured by its redshift Zs, and is given by 

^*-' ^ i^' <^' 

where dL{zs) is the luminosity distance and S{zs) is the K-corrected photon luminosity 
given by 

/•300keV 

S{zs)= / SE[{l + Zs)E]dE, (3) 

J50keV 

that is, the photon emission rate corrected for the redshifted bandpass. Evaluating this 
integral, we flnd 



^t ._ J Lo4.3x 10^(1 + zs 



2 - {1 + zs)-^ - 1(1 + z, 

D 



z. < 5 (4^ 



Lo4.3x 10''(l + 2s)"^5, z^>5 

where Lq is measured in ergs~^. Bursts with fluxes that are 5.5cr above background 
in either the 64 ms, 256 ms, or 1024 ms time bins trigger BATSE. The 1024 ms trigger 
threshold is most sensitive, so we adopt its flux limit, (/)min = 0.3 photon s~^ cm~^ for > 50% 
detection probability (Fishman et al. 1992). The flux limit (/>min deflnes the maximum 
redshift z^iax out to which bursts can be seen. 

The rate at which we observe GRBs that originate in redshift interval [^s, -2s + dz^] is 
given by 

dN no dV , , 

dzs {1 + Zs)dzs' 

The factor 1 + z^ comes from the gravitational time dilation of the burst rate, and dV/dzs 
is the comoving volume element. Using equations (2) and (5), we can calculate how 
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many bursts are observed with flux (f). The theoreticaUy predicted (F/Fmax) for this flux 
distribution is given by 

\ _ / I 1 \ _ _ (6) 



^max/ \ V^min/ / Jq'^''"'^ (dN/dZs) dZg 



The only missing ingredient required to carry out this computation is the specific 
choice of a cosmological model for the comoving volume element and the luminosity dis- 
tance dL- We will consider three Friedmann- Robertson- Walker cosmologies described by 
density and cosmological constant parameters (n,A) given by (1.0,0), (0.1,0), (0.1,0.9). 
Cosmologies with cosmological constant much more than 0.9 probably already are ex- 
cluded by the quasar lensing rate (Fukugita et al. 1992). The proper, luminosity, and 
angular distances are given by (c/. Carrol, Press, & Turner 1992) 

dmizs) = T7 

-"0 Jo 



{l + z)^{l + nz)-z{2 + z)A] ^'^ dz, 


(7a) 


dL = {l + Zs)dM, 


(7b) 


dA =dM/{l+Zs). 


(7c) 



In general, these distances must be computed numerically. The angular distance between 
two redshifts dA{zi, Zg) can be found by replacing the lower bound of integration in equation 
(7a) with the lower redshift. The comoving volume element is given by 

^ = ATTdl^dM = ^TTdl^ [{1 + Zs)H^ + fi^s) -Zs{2 + Zs)A] "'/' . (8) 

azs azg -no 

The quantitative effects of using an alternative, clumpy universe cosmology (Dyer & Roeder 
1973) are negligible for the moderate redshifts Zg ^ 1 out to which bursts can be seen 
(Fukugita et al 1992). 

The values of Lq, no, and -Zmax that reproduce the observed (V^/Knax) and burst rate 
are listed in Table 1. The parameters for O = 1, A = are in approximate agreement with 
those of Fenimore et al. (1993) (presumably the differences are a consequence of using a 
different spectrum). Our O = 0.1, A = 0.9 model also is consistent with that of Piran 
(1992), but no prior results are available for the O = 0.1, A = cosmology. We expect 
that a more careful analysis, for example including PVO data, would modify the results 
only slightly. We have not been more precise since burst modeling is not the aim of this 
paper, but is only one ingredient in the lensing calculation to follow. 

In Figure 1, the redshift distributions of bursts are shown. In fact, for each cosmology 
two curves are drawn. One curve is the computation of equation (5), and the second is the 
same computation, but showing the incompleteness as fluxes approach and fall below (j)uiin- 
We have found that the data in the BATSE efficiency table for 1024 ms triggers (Fishman 
et al. 1992) is fit very well by the function 

^w-{?-e-«-«/«, tit: («) 



GRB Lensing 5 

with (^0 = 0.21 photon s~^ cm~^ and a = 0.105 photon s~^ cm~^. We use equation (9) as 
the detection efficiency in this and aU subsequent calculations. The dotted vertical Unes 
are drawn at the values of -Zmax, where by definition the incompleteness is 50%. Note that 
some bursts originating at z^ > -Zmax will be seen when Poisson uncertainty in the photon 
rate brings the measured flux above (j)min- 

Band (1992, see also Petrosian 1993) has argued that failure to account for incom- 
pleteness near the flux threshold invalidates many comparisons between the observed and 
predicted {V/Vma.x)- We test the sensitivity of our parameters to incompleteness by includ- 
ing the detection efficiency function, equation (9) in our calculation of { .X ) . Multiplying 
the top and bottom integrands of equation 6by the detection efficiency and integrating 
the bursts between z = and z = -2max, we find (V^/Knax) = 0.27 instead of 0.32 for our 
O = 1, A = parameters. Furthermore, our and other models {e.g., TM93) reproduce 
the observed (V/Vmax) at thresholds well above the regime where incompleteness is rele- 
vant, but where cosmo logical effects are still important. Whether or not we include this 
incompleteness should modify the details of our cosmological model only slightly. 



3. THE RATE OF LENSING BY GALAXIES 

In this section we consider lensing by circular lens models. (See Blandford & Narayan 
1992, Narayan & Wallington 1992b for reviews of the physics of lensing.) Such lenses 
alway produce at least one image of magnification Mi greater than unity. This image is 
usually the brightest and most distant from the center of the lens. If there are multiple 
images, a second image with magnification M2 appears on the opposite side of the lens. 
This image can be very bright or very faint, depending on the lens model and the position 
of the source with respect to the lens. If the lens has a non-singular core, there can be a 
third image of magnification M3 in the core. This image is almost always very faint. 

The differential lensing rate of GRBs is given by 

dNi dN atotizs) 



dz^ dz^ 4-71" 



(10) 



^s 



where a"tot(-2s) is the total cross section of all lenses in the sky to produce an observable 
multiple burst from a source in the interval [^s, z^ + dzs]- If galaxies are the lensing masses, 
we can write the total cross section as 

(7tot(^s) = / dzi-— dL<5{L)a{zi,Zs,L). (11) 

Jo dz\ Jq 

The cross section for a single galaxy, a{zi, Zg, L), depends on the lens mass, which is directly 
related to lens luminosity. Weighting these cross sections by the galaxy luminosity function 
$(L), the inner integral sums over all lens masses at a single lens redshift zi, and the outer 
integral sums over all the lensing galaxies at each lens redshift. The lens population is 
characterized by the Schecter luminosity function, 

$(L)rfL = $*e-^/^* f ^"j ^, (12) 
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with parameters $* = 1.56 x 10~^/ifoo Mpc""^, [3 = —1.1, as assumed in Fukugita and 
Turner (1991, hereafter FT91). The dimensionless Hubble constant is defined by /iioo = 
ifo/lOOkms-^Mpc-^ 

The lensing cross section of a single galaxy of luminosity L to produce an observable 
multiple burst is 

a{zuz,,L)= [ dr,2nrse[M2{rs)(l){zs)]. (13) 

The variable Ts is the angular distance of a source to the lens axis. We have assumed we 
can relate the luminosity of a lens to its lens parameters. The value of Td-it, the outermost 
source radius that gives multiple images, depends on the particular lens model. Whether or 
not multiple images will be seen depends on the magnification M2 of the fainter image with 
photon fiux M2(j){zs)- If the detection efficiency €[M2(j){zs)] is below unity, the effective cross 
section is reduced below its maximum value vrr^^j^. (We note that Mao, 1992b assumed 
that the full isothermal cross section contributed. Blaes & Webster, 1992 did use the right 
selection criterion in their computation for point lenses.) Depending on the source position 
Ts, M2 can be greater or less than unity. If M2 is greater than unity, it is possible to see 
bursts beyond -2max7 but if it is less than unity, bursts may have to be considerably closer 
than -Zniax for two images to be detectable. 

In addition to the rate at which GRBs are lensed, we would like to know the charac- 
teristic time interval between the arrival of two images. The distribution of the time delay 
At between two bursts, which was computed first by Mao (1992b), is given by the formula 

dNi _ r,dN r^dV r^^^,^,da{z,,z,,L) 



dAt 



f^ , dN r= , dV [^ ^^ ^,^,da(zi,Zs,L) ,,,,,,,,, 



For any particular lens model, the time delay At can be related to the source position Ts, 
so that the differential cross section is 



da 2Txr^dr^ 



dAt ^drs 



(15) 



3.1. The Singular Isothermal Sphere (SIS) 

The singular isothermal lens (SIS) produces two images of the source whenever the source 
position lies inside the Einstein radius, Ts < r^, given by 

v^dA{zi,z^) 

where v is the line-of-sight velocity dispersion of the lensing galaxy. It is related to the 
lens luminosity via the Faber-Jackson or TuUy-Fisher relations, 

L^ \v^ J 
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Following FT91, the lens population is assumed to be comprised of E, SO, and S galaxies in 
the proportion 12%, 19%, and 69%, with corresponding parameters v^ = 225A/3/2kms~^, 
n = 4 for Es, v^, = 206 a/3/2 km s~^, n = 4 for SOs, and t)* = 144kms~^, n = 2.6 for 
Ss. The factor of a/3/2 in the velocity dispersion of Es and SOs accounts for isothermal 
dark halos surrounding more compact mass distributions (Turner, Ostriker, & Gott 1984; 
FT91). 

In the SIS lens, the magnification M2 is related to the source position by 

M2 = te/ts - 1. (18) 

For Ts > rE/2 (three quarters of the maximum cross section rcr'^), M2 is less than unity, 
and M2 goes to zero as Ts approaches ve, making the second image unobservable. In no 
case does a galaxy contribute the maximum cross section nr'^. 

Combining equations (10)-(13), (16)-(18) and evaluating the integrals numerically, 
we compute the lens rate as a function of source redshift z^ for each cosmology. The 
results are shown in Figure 2. For "standard" cosmological values O = 1, A = 0, we get a 
total multiple burst rate of 0.09 yr~^. This is quite a bit smaller than the value of 1.1 yr~^ 
predicted by Mao (1992b), both because he adopted a GRB model in which unlensed bursts 
could be seen out to redshifts -Zmax ~ 1-5 and because he used the maximum isothermal 
cross section. As expected (c/. Carrol, Press, & Turner 1992; Turner 1990; Fukugita, 
Futamase, & Kasai 1990), the finite A model gives the largest lensing rate, about 7 times 
larger than the standard model. 

The amplification bias is the ratio of the number of lensed bursts to the lensing rate 
computed as if bursts were not magnified. That is, it is the increase in the number of 
lensing events because sources intrinsically fainter than the flux threshold (^min are counted. 
We compute the hypothetical, no magnification lensing rate by replacing the detection 
frequency in equation (13) with e[0(2s)], i-e., assuming the magnification is unity. The 
cumulative burst rates and the amplification biases computed for the various cosmologies 
are presented in Table 2. The amplification biases are smaller than for quasar lensing. Two 
factors contribute to this. One is that the observed GRB flux distribution is not as steep, 
so that the increase in the number of observable sources resulting from amplification bias 
is not as large. The other, equally important factor is that amplification bias is computed 
using the magnification probability of the fainter second image. This differs from most 
computations of quasar lensing, where the magnification probability is computed either 
for the brighter image or for the total magnification {e.g., Turner, Ostriker & Gott 1984; 
FT91; but see Sasaki & Takahara 1993 for an exception). In either case, the magnifications 
of the quasars are always greater than unity, so that there is no loss of cross section for 
sources at ^ < -Zmax- For GRBs, however, we find a bias S ~ 1 for the fl = 1, A = model. 
Even though lensed GRBs can be seen to a somewhat higher redshift than unlensed bursts 
(compare Figs. 1 and 2), the loss of lensing cross section for small M2 almost perfectly 
offsets this increase. That is, even for sources at low redshifts, the second image will be 
too faint to see for some source positions that yield double images. In fact, biases can be 
less than unity (Sasaki & Takahara 1993), as in the O = 0.1, A = model. 

Even if our models of GRBs are more or less correct, reproducing the right Zmax, there 
are uncertainties in the lens rate resulting from uncertainties in the lens parameters. In 
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particular, since the lens rate varies as A^^ ~ (^*v^ ~ $*i'^(L/L*)^/"^, the 10% uncertainty 
in L^, 22% uncertainty in $*, and 10% uncertainty in i^n. {e..g., FT91; Mao 1992a) give 
nearly a 50% uncertainty in the computed lensing rates. Using a different observational 
luminosity function, Wallington and Narayan (1993) find a lensing rate about 40% larger 
than that given by the FT91 parameters. 

We next compute the statistical distribution of time delays. The time delay of a pair 
of images in a lens of Einstein radius te is given by 

At = f{z\, Zs, rE)rs, r^ < te (19) 

where we define the function 

./ N on , ^ ^a(0,^i) rfA(0,^s) 1^^. 

f{zi,z^,rE)=2{l + zi) —- -TE- (20) 

C dA[Z\,Zs) 



Then the differential cross section (eq. 15) is 

da 27rAt 



dM f 



(21) 



The corresponding source position is given by rs(At) = At/ f, and hence the magnification 
M2 of the second image can be calculated by equation (18). 

Substituting equation (21) into equation (14), we compute the distribution of time 
delays. The results are shown in Figure 3. The curves are divergent at small time delay 
due to the infinite number of low luminosity sources predicted by the luminosity function 
(eq. 12), but the total lensing rate, i.e., the integral under the curves, is finite and gives, 
of course, precisely the same cumulative lensing rate as the calculation for Figure 2. The 
median time delay is about 0.02 /i^qq yr for the two A = cosmologies and about 0.04 /i^qq yr 
for the O = 0.1, A = 0.9 cosmology. 

3.2. The Non-Singular Isothermal Sphere (NSIS) 

Realistic galaxy lenses probably have finite core radii. As is well-known, circular lenses 
with finite core radii can, in addition to the pair of images of a singular lens, produce a 
third image near the core of the lens. We adopt the usual non-singular isothermal lens 
model (NSIS, e.g., Hinshaw & Krauss 1987), for which we can calculate image positions 
and magnifications for any source position. The only difference between this calculation 
and the previous one is that we have replaced equation (18) with the appropriate relation 
for an NSIS lens. 

In previous investigations of quasar lensing (FT91; Mao 1992a), the importance of the 
core radius in galaxy lenses was based upon the observations of 42 elliptical galaxies by 
Lauer (1985). Lauer observed that 14 ellipticals had resolved cores, with half light radii 
that approximately follow the relation 

'■^'^^f^y\ (22) 



^1/2 ^ ^ 
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where r*/2(iA(-2i) = 0.16 kpc. The core radius, Tc, of the lens model is related to the half 

light radius by ri/2 = V^Tc- The remaining ellipticals had unresolved cores, so that we 
may regard the above values of ri/2 as upper limits. We assume all types of galaxies (E, 
SO, S) have core radii following the scaling of equation (22). The resulting lensing rates 
are reduced by about a factor of two, and this may be an overestimate if galaxy cores are 
often smaller than the above prediction. The lensing rates and amplification biases for the 
NSIS model are listed in Table 2. 

To compute how often the third core image is seen, we simply replace the magnification 
M2 in the detection efficiency computation with the magnification of the core image M3. 
We find that the third image will be observable in fewer than 1% of the double bursts. 
The third image will be visible only in the very rare, very luminous lenses where the core 
radius is large and the central image is only moderately demagnified. 

3.3. A Finite Mass Lens Model (QDEV) 

If galaxies, particularly ellipticals, do not have massive isothermal halos, lensing computa- 
tions based on a density profile containing finite mass may be more appropriate. For this 
analysis we adopt the quasi-de Vaucouleurs (QDEV) profile of Grossman and Saha (1994). 
This lens model is characterized by the defiection law 



AGM 

dA{Q,z\)c^ (re + rf 



«W = , ,. .^.2/^ , ^^2 • (23) 



This deflection law closely resembles the true de Vaucouleurs deflection law (Sanitt 1976) 
if the core radius Tc = 0.4re, where re is the half mass radius, but it is simpler to use in 
numerical calculations. 

We use the relations of Maoz and Rix (1993) to relate the luminosity of a lens to the 
lens parameters, core radius and mass. These relations are 

re = 0.4re = 0.4r: (^\ (24) 

and 

M fMV ( L 



where the parameters are r* = 4kpc, (M/L)* = 10, a = 1.2, b = 0.25. This lens can give 
three images, and as with the NSIS model, we are able to compute their positions and 
magnifications for any source position. Using the magnification of the second image M2 
in the calculation of detection efficiency, we find the lensing rate of this model. 

The lensing rate of GRBs, assuming all types of galaxy lenses can be modeled in this 
way, is given in Table 2. These rates are generally about half the SIS rate. The third core 
image is bright enough to be seen only about 0.1% of the time. The amplification biases 
are quite small. This is essentially because the cross section for multiple imaging is much 
larger than for the SIS lens (approaching the formally infinite cross section for a point 
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lens), but the cross section for two images to be bright enough to observe is similar to the 
SIS cross section. Although lensing permits sources somewhat fainter and more distant 
than -2max to be observed, the loss of lensing cross section associated with small values of 
M2 more than offsets this increase. 

The lensing rate of this model scales as Ni ~ ^*{M/L)L ~ ^*{M/L)* /L''^. The 22%, 
20%, and 10% uncertainties in $*, (M/L)*, and L* (FT91; Maoz & Rix 1993) combine to 
give about a 30% uncertainty in the computed lens rate. 

3.4. The BATSE Detection Rate 

The most pessimistic rate of observable double bursts in Table 2 is one every 25 years, 
while the most optimistic is one every 1.5 years. A standard O = 1, A = cosmology gives 
one every 10 years. Since GRO is likely to operate for at least several years, the odds are 
significant that there will be a lensed burst during this period if GRBs are cosmological. 

Unfortunately, BATSE is inefficient at detecting bursts. Due to Earth blockage re- 
sulting from CGRO's low orbit (450 km, Gehrels et al. 1993), BATSE only "sees" about 
2.67r steradians of the sky at one time, and the detectors have a duty cycle of only about 
55% (Fenimore 1993), giving a probability of 34% (Fishman et al. 1992) to see any single 
burst. Since the time delay between lensed bursts is typically much longer than the orbital 
period of the satellite (many days compared to ~100 minutes), the position of the satellite 
in orbit is essentially random with respect to burst arrival times. Thus, the odds of seeing 
both components of a double burst are 34% x 34% = 12%. In view of this instrumental 
inefficiency, if lenses are of the simple sort described above, it is very unlikely that multiple 
bursts will be seen. This instrumental inefficiency generally has been ignored in previous 
GRB lensing papers (an exception being Mao, 1992b), resulting in undue optimism {e.g., 
Blaes 1994) about the possibility of observing a lensing event in the lifetime of BATSE. 



4. LENSING BY AN ELLIPTICAL POTENTIAL 

About half of known multiply-imaged quasars show four, rather than two, bright images 
{e.g., Blandford & Narayan 1992). The simplest lens that can have such an image configu- 
ration is one with an elliptical gravitational potential. (See Narayan & Grossman, 1989 for 
a review of caustics, critical lines, and images in an elliptical lens.) If the requirement to 
observe a lensing event is for BATSE to detect two or more images of a quadruple burst, 
the instrumental detection efficiency may be much higher than the 12% of a double burst. 
In general, the probability of detecting two or more images out of n that are bright enough 
to be seen is given by the binomial probability 

P{>2\n) = J2{'i)f{'^-Pr-\ (26) 

i=2 ^ ^ 

where p = 34%. If only two images are bright enough to be seen, this gives P{> 2|2) = 12%; 
if three, P(> 2|3) = 27%; and if four, P(> 2|4) = 42%. Thus, if quadruple images are 
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sufficiently common, tlien tfie cfiance of detecting a lensing event may be better than 
predicted for circular lens models. 

To investigate the importance of elliptical lenses in GRB lensing, we model galaxy 
lenses by elliptical generalizations of the SIS potential, as used in the statistical lensing of 
quasars (Wallington & Narayan 1993). The scaled two-dimensional potential t(j (defined 
such that the deflection angle is ct = — V^) is given by 

i;{x, y) = TE [(1 + e)x^ + (1 - e)y^] '/' , (27) 

where e measures the lens ellipticity. We derive the statistics of elliptical lenses using 
a Monte Carlo simulation, where we randomly distribute isothermal lenses between the 
observer and source at z^ according to their luminosity distribution given by equation (12) 
and spatial distribution given by the co moving volume element. At each source redshift 
^s at 0.2 intervals, we randomly place 10,000 sources within a circle of radius 1.3 times 
the Einstein radius (since we are interested only in potentially observable multiple bursts) . 
The resulting statistical contribution to the time delay distribution of a pair of lensed 
bursts has a weight proportional to the area of the Einstein ring. 

All the images of a single source are located using an efficient grid search method, 
where we identify rough image locations on a coarse grid and iteratively subdivide to higher 
resolution only those grid cells which may contain an image. Each iteration increases the 
resolution by a factor of 10, and it goes 8 levels deep. Sometimes several neighboring pixels 
are marked as containing an image, so as a final step, we compute the x and y components 
of the lens equation {e.g., Xs = x — a{x, y)) at the four corners of each pixel and keep only 
pixels where the source position is bounded by the four corner values. We use the four 
corner values of Xg and ys to interpolate x and y to a refined image position. 

The results for a lens of ellipticity e = 0.2 are shown in Figure 4, where we plot as a 
bold line the ratio of the cross section of the elliptical lens for observable multiple imaging, 
a"e(2-|-), to the cross section of the circular lens for double imaging, cro(2). The cross section 
a"e(2+) is divided into the contributions of two-image bursts, three-image bursts, and four- 
image bursts, i.e., cre(2-|-) = a ^{2) -f cre(3) -|- (7^[A). Analyzing the data of Huchra et al. 
(1983), Wallington and Narayan (1993) show that few galaxies have ellipticities exceeding 
e = 0.3. They find a mean ellipticity somewhat less than e = 0.2, and thus, using e = 0.2 
in Figure 4 is roughly representative of a typical galaxy lens. 

The caustic curves of an elliptical lens divide the source plane into regions capable of 
producing one, three, and five images. The five image region is bounded by a diamond- 
shaped astroid, and the associated cross section for producing five images is a^. The 
highest magnification events are associated with sources near folds or cusps on the caustic. 
The three image region is entirely outside this astroid (as long as the core radius is small 
enough), and is bounded on the outside by an ellipse with mean radius te- The associated 
cross section is (J3, and bright images associated with appear on opposite sides of the lens. 
For lenses with small core radius, the central image of a three or five image lens is highly 
demagnified, and since we use a singular lens, the central image is always absent. Thus, 
the as and sigma^ cross sections produce two and four images, respectively. 
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The elliptical lens produces fewer multiple imaging events with images above the 
detector threshold than does the circular lens. This is essentially because, unlike a circular 
lens which produces high magnification images in pairs, an elliptical lens can produce one 
(or more) bright image of a source near a cusp, while leaving its counter-image magnified 
much less {e.g., Narayan & Grossman 1989). Highly magnified pairs or triplets of images 
are likely to be mergers across folds and cusps. Counter-images with lower magnification 
will not be seen at redshifts far beyond -Zmax, and indeed there are essentially no quadruple 
images in Figure 4 beyond ^s ~ 1-2. 

For a lens with e = 0.2 and small core radius such that the diamond caustic (defining 
the five-image region) lies entirely within the the radial caustic (bounding the three-image 
region), the three-image cross section a^ is greater than the five-image cross section a^ 
at low and moderate magnifications (Wallington & Narayan 1993). Images and counter- 
images can be seen if a source is at low redshift, and, therefore, pairs of images dominate the 
multiple image cross section a{2+). The three-image cross section varies as a^ ~ M~^'^, 
while the five-image cross section varies as a^ ~ M~^, so that at large magnifications, the 
five-image cross section dominates over the three-image cross section. The transition to a^ 
dominance is at magnification M ~ 10 (Wallington & Narayan 1993). Thus, three- and 
four-image bursts become increasingly more important at larger redshifts, where bursts 
must be magnified more just to be seen. (We associate both three- and four-image bursts 
with a^ since they arise mainly from cusp imaging. Two-image bursts are mainly associated 
with (73 when they are at low redshift. At high redshift, however, the double images are 
mergers across folds, rather than image/counter-image pairs, so that these pairs are also 
associated with a 5.) However, due to the low magnification bias, very high magnification 
events are not very common, and three- and four-image bursts are never seen to dominate 
over two-image bursts, unlike quasar imaging. 

Although the absolute rate of observable multiple imaging by an elliptical potential 
is reduced by factor a^{2+)/ao{2), it is possible that the rate at which lensing events are 
detected is increased. The enhancement of the rate at which lensing events are detected is 
given by the multiple image cross sections, weighted by the BATSE detection probability, 
according to 

iV, ^ 0.12a,(2)+ 0.27(7,(3) + 0.42(7,(4) 

A^o 0.12(70(2) • ^ ' 

This ratio is drawn as the dashed line in Figure 4. At the low and intermediate redshifts 
that give quadruple bursts, the detection rate is enhanced above the circular lens rate. 
For the most probable source redshift of Zs ~ 0.8 (see Fig. 2, O = 1, A = curve), the 
detection rate is about 30% larger than for a circular lens. If e = 0.1, this enhancement 
is about 15%, and if e = 0.3, it is about 70%. Averaged over all source redshifts, the 
enhancement of the lensing rate clearly is less. Thus, BATSE' s efficiency to detect lensing 
events is only slightly better than 12%, being perhaps around 15%. 

In Figure 5 we show the time delay distributions for the two Monte Carlo simulations 
presented in Figure 4. The e = distribution closely resembles the O = 1, A = model in 
Figure 3, except for the absence of a divergence at small At, resulting from the absence of 
lenses with L < O.IL^, in the simulation. The e = 0.2 distribution shows many more short 
time delay bursts, resulting from images merging across folds. (This distribution includes 
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all visible pairs in a multiple image burst. For example, a quadruple burst has six time 
delays.) The median time delay in this distribution is about three times shorter than in 
the e = distribution. (Both medians here are systematically too large due to the absence 
of low mass lenses, but if we calibrate the e = median with the 0.02 /i^qq year interval 
computed from the analytic model, we can get an idea of the true median of the e = 0.2 
distribution.) Thus, the time delays associated with triple and quadruple imaging may 
often be less than a day, and most often will be associated with the faintest bursts which 
are seen only at high magnification. 



5. DISCUSSION AND CONCLUSION 

We have presented three new cosmological models of the GRB population, derived using a 
double power-law model of burst spectra (TM93), and have computed the rate of double 
imaging of GRBs by the known population of galaxies, using three different models of 
galaxy lenses. The highest lensing rate is derived assuming a large cosmological constant 
and a singular isothermal lens, while the lowest rate has zero cosmological constant and 
a lens with finite core radius. We predict there should be one lensing event above the 
BATSE detection threshold every 1.5-25 years. The median time delay between the pair 
of images is ~ 0.02 years. Unfortunately, the instrumental inefficiency associated with 
BATSE's duty cycle (12% probability of detecting both components of a double burst) is 
so severe that seeing a lensing event from a galaxy lens is highly unlikely. 

The instrumental inefficiency of observing a lensing event is improved if four com- 
ponents of a burst, such as produced by an elliptical potential, are above the BATSE 
threshold. In this case there is a 42% chance of detecting at least two of the four images. 
Monte Carlo simulations of lensing by elliptical potentials show, however, that quadruple 
imaging is less common than for quasar lensing, essentially because the amplification bias 
of lensed bursts is smaller. Consequently, the rate at which lensed bursts will be detected 
is increased by at most about 20%. 

Despite our pessimistic prediction for the detection rate of lensing events, we believe 
it is well worth the effort to study the BATSE light curves for possible lensing events. 
A single lensing event, albeit improbable, would verify the cosmological origin of GRBs. 
Furthermore, even a single detection would suggest one of two exciting possibilities. It 
might indicate that the universe has a nonzero cosmological constant, giving a larger 
lensing rate than the "standard" (O = 1, A = 0) model. An alternative explanation 
could be the existence of a significant cosmological density of compact massive lenses [e.g., 
10^ Mq black holes), in which case the lensing rate could be much higher than for galaxy 
lensing (Mao 1992b; Blaes & Webster 1992). The absence of lensed bursts can be used to 
set constraints on the density of such massive compact objects (Nemiroff et al. 1993). 

Nowak and Grossman (1994) have proposed statistical tests for comparing noisy light 
curves and find that many bursts of comparable duration, but distinct origin can be dis- 
tinguished by their methods. Their methods are extremely reliable for the 41% of bursts 
with 4'/4>uiin > 3, but for fainter bursts noise often will be sufficient to obscure the differ- 
ences in the light curves. Thus, consistent pairs of light curves among the brighter bursts 
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would be good lensing candidates. Of course, one can only demonstrate that two bursts 
are different, but not that they have the same source. The frequency of false positives, that 
is, error measures that are consistent with the lensing hypothesis for two bursts that are 
actually different, could be investigated by applying the statistical tests to bursts which 
are known to be distinct {e.g., because they appeared in different parts of the sky). If false 
positives are very rare, then a pair of consistent light curves may constitute strong evidence 
of lensing. Calibration of the false positive rate has not been undertaken, however, due to 
the enormous number of comparisons (thousands) involved. Nowak and Grossman have 
examined 260 bursts preceding March 5, 1992 in the BATSE catalogue (Fishman et al. 
1992), and identified 27 pairs of bursts of comparable duration whose positional error boxes 
overlap. As expected, there are no candidate lensing events among the bright {4>/4>iQin) 
bursts, and only one fainter statistically consistent pair. However, the duration of this pair 
of bursts is so short that there are not enough time bins in the light curves for high quality 
statistical comparison. Thus, to date there are no good lensing candidates in the data. 

We emphasize that the improbability of detecting a lensing event is not physical, but 
mainly instrumental. An experiment designed with lensing in mind could greatly increase 
the odds of a successful lensing detection. Most obviously, a satellite in a higher orbit 
could have a duty cycle near 100%, both because of smaller Earth blockage and because 
the detectors could operate continuously. Then, instead of 12% instrumental efficiency 
to detect both components of a lensed burst, there would be nearly 100% instrumental 
efficiency. We also note that three BATSE detectors in a comparable low orbit could also 
give 100% efficiency (almost a ten-fold increase in efficiency for three times the cost!). 

TM93 have proposed that a "BATSE detector" with 18 times greater sensitivity could 
detect enough distant bursts in three years to distinguish between go = 0.1 and qo = 0.5 
cosmologies. In this case the lensing rate would be about 4 times larger, giving one 
lensed burst every 0.4-6 years, with 2.7 years being the mean for the O = 1, A = 
cosmology. If GRBs are cosmological in origin, an appropriately designed experiment with 
high instrumental efficiency would probably detect lensing events. Failure to detect any 
candidate lensing events within 4.6 times the mean burst interval (say 4.6 x 2.7 ~ 12 years) 
would exclude the cosmological scenario at 99% confidence. 



The authors are indebted to Man Hoi Lee and Chris Thompson for discussions and 
for pointing out certain references during the course of this project. 
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TABLE 1 
Cosmological GRB Models 



Q A Lo/ifoo^rgs ^ noh^oojr ^ Mpc" 



1.0 


0.0 


1.1 X 10^° 


7.8 X 10-s 


0.83 


0.1 


0.0 


2.0 X 10^° 


4.2 X 10-s 


0.94 


0.1 


0.9 


6.0 X 10^0 


1.3 X 10-s 


1.13 



TABLE 2 
Cumulative Lensing Rates and Amplification Biases 

SIS NSIS QDEV 

n A Ni{yr-^) B Ni{yr-^) B Ni{yr-^) B 



1.0 


0.0 


0.09 


1.02 


0.04 


1.91 


0.04 


0.01 


0.1 


0.0 


0.17 


0.98 


0.08 


0.48 


0.07 


0.0005 


0.1 


0.9 


0.67 


9.94 


0.35 


1.56 


0.20 


0.0004 
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FIGURE CAPTIONS 

1. The redshift distribution of GRB as computed from equation (5), with and without 
correction for incompleteness at faint fluxes near the detection threshold. The curves 
are labeled by the cosmological parameters (n,A). The dotted vertical lines indicate the 
values of ^max, which by definition, correspond to 50% completeness. The integrals under 
the curves allowing for incompleteness give 800 bursts per year. 

2. The differential lensing rate as a function of source redshift ^s for the SIS lens model. 
The curves are labeled by the cosmological parameters (O, A), and the cumulative lensing 
rates A^; [i.e., the integrals under the three curves) are given at upper right. 

3. The lensing rate as a function of time delay At for the SIS lens model. The curves are 
labeled by the cosmological parameters (O, A), and the cumulative lensing rates A^^ (i.e., 
the integrals under the three curves) are given at upper right. The median time delays are 
about 0.02 /i^QQ years (7 days). 

4. The rate of multiple imaging by an elliptical potential (e = 0.2), measured by cross 
section a^{2+), compared to the rate of double imaging by a circular potential, measured 
by cross section cro(2), is drawn as the solid, bold line. The multiple imaging cross sec- 
tion is divided into the contributions of double, triple, and quadruple images. (That is, 
the bottom curve is cre(2), the next one is cre(2) + cre(3), and the next one, cre(2-t-), is 
ae(2) + ae(3) + a^{4).) Double images dominate at low redshift, but triple and quadruple 
images are increasingly important at higher redshift. Due to the low amplification bias, 
they never dominate, however. Accounting for BATSE's detection efficiency, we show the 
enhancement of the detection rate, N^/Nq, as the dashed line. The elliptical potential gives 
an increased detection rate of at most a few tens of per cent. At redshifts Zg ^ 2.5, fewer 
than 10/10000 sources in the Monte Carlo simulation have observable multiple images, 
and numerical noise dominates the results. 

5. This distribution of time delays between pairs of images derived from the Monte Carlo 
simulation for e = (solid histogram) and e = 0.2 (dashed histogram). The curves should 
be compared to the O = 1, A = model in Fig. 3. The curves do not diverge at small 
At because of the absence of lensing galaxies fainter than L = O.IL* in the simulation. 
Images merging across folds produce many more short time delays in the e = 0.2 model 
than in the e = model, giving a median time delay about three times shorter. 
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Figure 1 
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Figure 2 
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Figure 3 
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Figure 4 
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Figure 5 
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